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ABSTRACT 



The Fast Fourier Transform is an algorithm for the computation 
of Discrete Fourier Transforms in less time than allowed by any other 
algorithm available. The use of special purpose digital machines 
to reduce those times even further is of interest for real time 
spectral analysis. The main principles of Fast Fourier Transforms 
are presented. The design of a full-parallel eight sample processor 
is presented as a point of reference for comparison with serial and 
serial-parallel hybrid machines. Carry-Save Addition is introduced 
and used as the primary arithmetic logic. 
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th 

The k sample of a sampled time signal 



7 










fi 











1 • Introduction 



In presenting a procedure that will determine the complex Fourier 

coefficients of a sampled signal, it is first necessary to introduce 

the discrete Fourier transform (DFT) and its inverse. The DFT is 

defined by 
N-1 

F(n) = ^ Z(k)exp(-j2nnk/N) k = 0,l,2, N-1 (1) 

k=0 

where F(n) is the n^^ coefficient of the DFT and Z(k) the k^^ sample 
of a N sample signal. It is assumed that the N samples are equally 
spaced and taken at a frequency that is at least twice the highest 
frequency component of a bandlimited signal. Using the principles of 
orthogonality it can be shown that the inverse of the DFT is 

N-1 

Z(k) = Y F(n)exp( j2TTnk/N) n = 0,1,2, N-1 (2) 

k=0 

Since the procedure uses successive reductions of a finite series* 
length by a factor of two, the number of samples, N, must be an integer 
power, p, of 2 giving N = 2^. It is evident at this point that p 
reductions or passes will be necessary to obtain a non-serial solution 
for a given coefficient. 

Partitioning equation (1) at the halfway point, it appears as 
N/2-1 N-1 

F(n) = Y Z(k)exp(- j2nnk/N) + ^ Z (k)exp(- j2nnk/N) (3) 

k=0 k=N/2 

substituting k = m + N/2, in the second half of the equation yields 
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( 4 ) 



N/2-1 N/2-1 

F(n) = y Z (k)exp(- j2TTnk/N)+ ) Z(m+r)exp[- j2nn(m+r) /N] 

k=0 m=0 2 2 



and combining under a single summation 



N/2-1 

F(n)= ' Z (k)exp(- j2rrnk/N)+Z (k+r)exp[ - j2nn(k-)7-'' /"”) 1 

1^=0 2 , , 



(5) 



The common exponential multiplier may be factored and after reduction 
of the internal exponential one obtains 



N/2-1 

F(n)= [Z (k)+Z(k-W/2)exp(- jnn)]exp(- j2nnk/N) 

k=0 



( 6 ) 



where 



+1 n even 



exp(-jnn) = 



• 1 n odd 



giving 

N_/2-l 

F(n) = [Z(k)+Z(k-W/2)]exp(- j2TTnk/N) 

k=0 



n even 



(7) 



N/2-1 

F(n) = [Z(k)-Z(k4N/2)]exp(- j2nnk/N) n odd 

k=0 



( 8 ) 



Initiating a new pass, equations (7) and (8) are partitioned at the 
half-range point 



F(n) = 



N/4-1 

, [Z (k)+Z (k4U/2) ]exp(- j2rrnk/N) 
k=0 

N/2-1 

+ y [Z(k)+Z(kdN/2)]exp(- j2TTnk/N) n even 

k^/4 



(9) 
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F(n) 



N/4-1 

y [Z(k)-Z(k+N/2)]exp(- j2rrnk/N) 
k=0 

N/2-1 

+ y [Z(k)-Z(k+N/2)]exp(- j2TTnk/N) n odd 
lc=N/4 



( 10 ) 



substituting m = k - N/4, k = m + N/4, in the second summations 



N/4-1 

F(n) = Y [Z(k)+Z(k-W/2)]exp(- j2nnk/N) 
k=0 

N/4-1 

+ y [Z (m+N/4)+Z (m+3N/2)]exp[ - j2nn(m-fN/4) /N] n even (11) 



N/4-1 

F(n) = y [Z(k)-Z(k4N/2)]exp(- j2TTnk/N) 
k=0 

N/4-1 

+ j [Z(m+N/4)-Z(m+3N/2)]exp[- j2TTn(m+N/4)/N] n odd (12) 
m=0 



Combining under a single summation and extracting the common exponentials 
gives 



F(n) 



y j [Z(k)+Z(k-»I/2)]+[Z(k-W/4)+Z(k+3N/4)]exp(- jTTn/2) 

1^0 ' 



X exp(- j2nnk/N) n even (13) 

N/4-1 

F(n) = Y [Z(k)-Z(k-tN/2)]+[Z(k+N/4)-Z(k+3N/4)]exp(- jTTn/2) 
k=0 ' 



X exp(- j2TTnk/N) n odd 



(14) 



where 
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1 


nil 


even 


-1 


n/2 


odd 


exp(-jnn/2) = ( 






] -j 


(n-l)/2 


even 




(n-l)/2 


odd 



After two passes the coefficients are defined by 



F(n) 




f I [Z(k)+Z(k+N/2)]+[Z(k-W/4)+Z(k+3N/4)] ! exp(- j2nnk/N) 
1^0 ’ ■ 



n/2 even 



(15) 



N/4-1 . 

F(n) = y ; [Z(k)+Z(k-W/2)]-[Z(k-™/4)+Z(k+3N/4)] 'i exp(- j2nnk/N) 

1^0 



n/2 odd 



(16) 



N/4-1 , 

F(n) = r I [Z(k)-Z(k-W/2)]- j[Z(k-HN/4)-Z(k+3N/4)] . exp(- j2nnk/N) 
k=0 

(n-l)/2 even (17) 



" L i [Z(k)-Z(k-™/2)]+j[Z(k-ttJ/4)-Z(k+3N/4)] j exp(- j2nnk/N) 
k=0 

(n-l)/2 odd (18) 



General characteristics of each reduction are now apparent. They are 

1. For the pass a complex multiplier of the form exp (-j2rrn/2^) 

is generated. 

th 

2. The limits of the summation after the u pass are k=0 to 
k=(N/2^)-l. 
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3. For each reduction beginning with some combination of samples, 

T(k), the multiplicand of the complex multiplier is T(k‘W/2^). 

th 

4. The general form after the u pass is 
N/2^" ^ 

F(n)= y LT(k)-f«xp(- j2TTn/2“) T(k-W/2'^)]exp(- j2TTnk/N) 
k=0 

5. When u=p, N=2^, the limits of the summation are k=0 to k=0; 
the exponential exp(- j2TTnk/N) becomes unity and the non-serial solution 
of F(n) is explicit, 

r d 

Now, continuing the procedure for the 3 pass, u=3 : 

a. the exponential multiplier is 

3 

exp(-j2rrn/2 ) = exp(- jrrn/4) ; 

3 

b. the limits of summation are k=0 to k=(N/2 )-l=N/8-l; 

3 

c. the sample shift is N/2 = N/8; 

d. the general form after the third pass is 

N/8-1 

F(n)= y [T(k)+exp(- jrrn/4)T(kfN/8)]exp(- j2rrnk/N) . 

k=0 

e. if N=8, 

F(n) = T(0)+exp(- jTTn/4)T(N/8) . 

Recalling equations (15) through (18) the following substitutions 
are made to simply bookkeeping: 

A(k)=Z(k)+Z(k“W/2) B(k)=Z(k-W/4)+Z(k+3N/4) 

C(k)=Z(k)-Z(k+N/2) D(k)=Z(k-W/4)-Z(k+3N/4) 
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Also, in any set of equations having common summations, the limits of 
the summations will only be shown for the first equation of the set. 

Rewriting (15-18) yields 



N/4-1 

F(n) = y [A(k)+B(k)]exp(- j2TTnk/N) 

k=0 


n/2 


even 


(15a) 


F(n) = 5] [A(k)-B(k)]exp(- j2Trnk/N) 


nil 


odd 


(16a) 


F(n) = [C(k)-jD(k)]exp(-j2rrnk/N) 


(n-l)/2 


even 


(17a) 


F(n) = [C(k)+jD(k) ]exp(- j2TTnk/N) 


(n-l)/2 


odd 


(18a) 



After the third pass, the resulting equations are 



N/4-1 , 

F(n) = y [A(k)+B(k)]+ exp(-jTTn/4)[A(k+N/8)+B(k+N/8)] exp(-j2rmk/N) 

k=0 ‘ ' 

n/2 even (19) 



F(n) = [A(k)-B(k)]+ exp(- jrrn/4)[A(k+N/8)-B(k+N/8)] | exp(- j2Trnk/N) 

n/2 odd (20) 



F(n) = X I [C(k)- jD(k)]+ exp(- jTTn/4)[C(k-W/8)- jD(k-m/8)] | exp(- j2nnk/N) 

(n-l)/2 even (21) 



F(n) = [C(k)+jD(k)]+ exp(- jrrn/4)[C(k+N/8)+jD(k+N/8)] | exp(- j2Trnk/N) 

(n-l)/2 odd (22) 
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noting 



exp(- jTTn/4) 



+1 


n/4 


even 


-1 


n/4 


odd 


+.707(l-j) 


(n-l)/4 


even 


-.707(l-j) 


(n-l)/4 


odd 


- j 


(n-2)/4 


even 


+j 


(n-2)/4 


odd 


-.707(l+j) 


(n-3)/4 


even 


+.707 (1+j) 


(n-3)/4 


odd 



Performing the complex multiplication yields 



F(n) 



N/8-1 

T i [A(k)+B(k)]+[A(k-W/8)+B(k-H'I/8)] exp(- j2rrnk/N) 
k=0 ' ' 

n/4 even 



F(n) = Y. I [A(k)-fB(k)]-[A(k-W/8)+B(k-«^/8)] | exp(- j2Trnk/N) 

n/4 odd 



F(n) 



Y j C(k)+.707[C(k+N/8)-D(k-W/8)] I 



-j D(k)+.707[C(k-«J/8)-H)(k-«I/8)] j 



exp(- j2rrnk/N) 



(n- 1) /4 even 



F(n) = i C(k)-.707[C(k+N/8)-D(k-«l/8)] j 

-j : D(k)-.707[C(k-«J/8)-m(k-W/8)] Jexp(- j2TTnk/N) 

(n-l)/4 odd 

F(n) = Y j [A(k)-B(k)]- j[A(k+N/8)-B(k+N/8)] 1 exp(- j2TTnk/N) 

(n-2)/4 even 



(23) 



(24) 



(25) 



(26) 



(27) 
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( 28 ) 



F(n) = ^ i [A(k)-B(k)]+j[A(k+N/8)-B(k4N/8)] j exp(- j2TTnk/N) 



(n-2)/4 odd 



F(n) = j C(k)-,707[C(k-W/8)-D(k-W/8)l j 

+j j D(k)-.707[C(k-W/8)+D(k+N/8)] ) jexp(- j2TTnk/N) 

(n-3)/4 even (29) 

F(n) = 5:^, C(k)+.707[C(k-W/8)-D(k-W/8)] j 

+j j D(k)+.7 07[C(k+N/8)+D(k-W/8)] j | exp(- j2TTnk/N) 

-i 

(n-3)/4 odd (30) 

Substitution for the functions A,B,C, and D will produce explicit statements 
of the transform equations after three passes. 

Further arithmetic development is avoided for brevity with the state- 
ment of the transform equation sets for the 8 sample (N=8) and 16 sample 
(N=16) cases in Appendices I and II, respectively. Appendix III, in turn, 
is an illustrative problem in which the spectral lines of a signal are 
determined and subsequently the originating signal reproduced via the 
described transform and its inverse. 
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The procedure just presented is indeed a very simple example 
of the Fast Fourier Transform (FFT) , an algorithm that allows the 
computation of a time series DFT with fewer actual arithmetic 
operations than other algorithms available. An arithmetic operation 
is defined, here, as a complex multiplication followed by a complex 
addition. 

Noting that the straight forward method of computing the DFT 

2 

requires N operations, Cooley and Tukey [l] showed that less than 
2Nlog2N operations are required when the FFT algorithm is used. 

They also showed that this algorithm requires no more data storage 
than the storage required for the initial samples, assuming the 
initial samples are complex. The particular method demonstrated 
herein, however, was first shown by Sande [2] and later referred to 
as "Decimation in Frequency" by Cochran, et al, [3] because of the 
characteristic divisions of the F(n)*s after each pass. 

In Fig. 1 a flow diagram for the equations of Appendix I, the 
eight sample example, is shown. To be sure, this diagram is not 
entirely general; the assumption of real inputs avoids continual 
complex additions and allows absolute segregation of the real and 
imaginary parts of the computed coefficients. Such an assumption 
has many real world applications and succeeds in reducing the 
computation time and circuitry by at least a factor of two. 

A flow diagram for the computation of the even numbered coef- 
ficients of the 16 sample example (equations of Appendix II) is 
presented in Fig. 2. It is interesting to note that this is, in 
fact, the same data flow as shown in Fig. 1, with the exception 
of the first addition cycle. If the first subscript of each basic 
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SAMPLE FIRST SECOND THIRD 

SUBSCRIPT PASS PASS PASS 

Fin} 

REAL IMAG 
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Figure 1. The flow diagram for the coefficients of the eight sample problem. As in all 
following flow diagrams, the samples are referred to by their subscripts. The 
value next to an arrow on a flow line references a multiplication, while a node 
infers a sum. 
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SAMPLE FIRST SECOND THIRD F(n) 
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Figure 3. The flow diagram for the odd coefficientc of the sixteen sample problem. 
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Figure 4. The characteristic reduction of the Fast Fourier Transform by pre-multiplication 
of the basic couplets with appropriate factors before processing through the next 
lower order transform. 



couplet in the first pass column in Fig. 2 were used to define the 
couplet, the last statement becomes even more obvious. Figure 3 
completes the flow diagram for the 16 sample example by presenting 
the flow for the odd coefficients. Notably, this type of flow 
does not take full advantage of the FFT properties. 

The pre-multiplication of difference couplets by appropriate 
exponential factors and then continuing through the next lower order 
process was shown graphically in [3] for the eight sample case. 

Figure 4 extends that representation to the 16 sample case and 
illustrates the ease with which general purpose computers might 
handle such computations. Comparing Figures 3 and 4, the inherent 
trade-off between special purpose and general purpose implementations 
of a particular function is illustrated. The use of recursive proce- 
dures and the availability of complex arithmetic makes general purpose 
computation of the Fig. 4 process most advantageous. Whereas the 
elimination of complex multiplications and complex input data in the 
Fig. 3 process produces a more cumbersome set of equations; equations 
which are, however, more easily adapted to special purpose computation. 

It is the purpose of this work to investigate the applicabilities 
of special purpose computers to Fast Fourier Transforms, The main 
criterion of comparison between special purpose and general purpose 
computations will be time versus circuit complexity, where complexity 
includes the number of components, the size of memory, and the size 
of the control unit. 

A completely parallel, 8 "real" sample, 8 spectral line processor 
will be used as a vehicle for the "fast" special purpose machine. 
Extrapolations to greater numbers of sampled inputs and the subsequent 



22 



increase in spectral lines will be made with respect to the increased 
complexity and cost. Further, an investigation of hybrid sequential- 
parallel machines will be presented. 
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2 . Parallel Processing 



The data flow diagrams of the eight sample example is shown in 
Fig. 1 and the explicit set of equations given in Appendix I. 

As the word "parallel" has many connotations when referring 
to computers, "full parallel" computation implies the simultaneous 
computation of all equations of the set. This does not, however, 
mean that the computation of each equation in all passes of the 
FFT must be independent. In fact, independent computation, when 
applied to the particular set of equation being examined, would be 
ridiculous. The computation in each pass is done in parallel. 

Examining Fig. 1, six basic combinations of four samples (quartet) 
and two basic couplets are found at the end of the second pass. Fo- 
cusing attention on the quartets, for the present, the principle of 
Carry-Save Adders (CSA) is utilized. As shown in Appendix IV, CSA 
techniques may be used to obtain exceptionally fast addition of 
several arguments through logic rather than cyclic operations as in 
conventional adders. It might be added that CSA may also be used to 
implement extremely fast multiplication processes. Six 4-argument 
CSA*s will be used at the front of the purposed processor. 

At this point, the real parts of F(2) and F(6) and the imaginary 
part of F(6) are available. The addition of one more rank of CSA*s 
to the imaginary F(6) processor provides the imaginary part of F(2), 
the negative of the imaginary part of F(6). F(0) and F(4) are obtained 

by crossfeeding the results of the all even samples and all odd samples 
summers into two additional two-argument CSA's. The odd- sum must be 
complemented in one case, producing F(4); the non- complemented case 
yields F(0) . 
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igure . The block diagram of the full-parallel eight sample processor. The numerical 
signs in parentheses indicate the action within the CSA units. The absence of 
parentheses indicates a positive input. 



To obtain the odd coefficients, the results of the two remain- 
ing quartet CSA*s must be multiplied by the constant . .70711. Multi- 
plication is a shift and add process, dependent on the position of 
"ones" in the binary multiplier. Since .70711^^ = .lOllOlOlOO^ the 
shift would be to the right and five additions necessitated. 

For example, the multiplication of an eight bit quantity M by 
.70711 would be represented as: 



M, M, M 

76543210 

iiQ 



MMMMMMMM 

U7 ii^ ii^ ii^ ii^ ii^ ii^ 



MMMMMMMM 

7 6 5 4 3 2 1 0 



PPPPPPPPPPPPPPPP 
7 6 5 4 3 2 1 0 -1 -2 -3 -4 -5 -6 -7 -8 



A CSA multiplier for this example would promulgate the carrys from the 
-6 digit forward and require a maximum of five arguments (the P^ bit). 
The maximum number of delays for multiplication by the given constant 
would then be 5+k for a k bit multiplicand. 

Performing the multiplication in two Carry- Save Adder (CSA) 
multipliers then allows the crossfeeding directly into four 3-argument 
CSA adders. This gives the respective combination for the real parts 
of the odd coefficients and the imaginary parts of F(3) and F(7). 
Complementing the imaginary parts of F(3) and F(7) by an additional 
CSA rank on each of the sums yields the imaginary part of F(l) and 



F(5). 



Having outlined the complete processor verbally, the block diagram 
is presented in Fig. 5. Note, the computation of the imaginary parts 
of F(l) and F(5) are the most involved; therefore, the time analysis is 
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based on those computations. For the sake of real time estimation, 
a 12-bit word and the availability of the 20 nsec adders mentioned in 
Appendix IV are assumed. The number of delays to the final stage is 
taken from Appendix IV and 11 ripple delays are assumed for the 12- 
bit word. 

The initial four argument processor required three delays, the 
multiplication four delays, the final three argument addition two 
delays, and the complement one delay for a total of 10 delays to the 
final stage. Adding 11 ripple delays gives 21 units of delay for a 
420 nsec total time of computation for the eight spectral lines of 
a signal. 

A maximum of k - 1 full adders are required for each bit of a 
k argument CSA. A measure of complexity, the total number of adders, 
may be calculated. Thirty-six adders are required for each of the 
six, 4-argument, 12-bit units for a total of 246 adders for all six 
units. For each of the two argument summations and the three negation 
operations, 12 adders are required for a total of 60. For each of the 
four 3-argument units following the multipliers, 24 adders are required 
for a total of 96. Finally, for each multiplier 45 adders are needed 
giving a 90 adder total. The processor then requires a grand total 
of 492 adders. A completely parallel computation of the eight Fourier 
coefficients takes less than .5 usee, using five hundred 20 nsec 
full adders. 

Since the initial information may be stored in the same location 
as the eventual answers, 192 bits of immediate access memory is needed. 
No control units other than a .5 usee clock and a memory read and write 
logic are required. If J-K flip-flop memory is used, only one AND 
gate is needed for control of the memory. 
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Extension of the methods of computation shown for the eight 
sample example to a sixteen sample example may be accomplished with 
one basic assumption; an average of five additions is required for 
each constant multiplier. 

Referring to Fig. 2, the simplest method of obtaining the even 
coefficients would be to provide eight parallel adders in front of 
an eight sample processor. For the specified 12-bit word, this 
requires an added 96 adders to the basic 492 and would add 20 nsec 
to the computation time. 

Computation of the odd coefficients becomes slightly involved. 

It is found that 1076 adders are required and an approximate compu- 
tation time of 500 nsecs achieved (based on the multiplication 
assumption stated above). 

Thus a 16-sample processor would require a total of 1664 adders 
and would produce the 16 complex spectral lines in 500 nsec. It is 
quite evident that extension of the **full parallel" computation to 
higher number of imputs becomes extremely expensive in terms of the 
number of adders required. However, the computation times are 
extremely fast. It is estimated that almost two million adders would 
be required for the 1024-sample case. The computations would still 
be accomplished in less than 2 usee. 

It must be remembered that the figures given in this section 
are for real input values, not complex. To further extend the 
processors to handle complex input signals would require an approximate 
257o increase in adders and a 107o increase in processing time. Thus, 
the eight and sixteen sample processors would require approximately 
600 and 2000 adders respectively for the computation of coefficients 
in .5 usee. 
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3. Hybrid Processing 



A FORTRAN IV program used to compute the Fast Fourier Transforms 
on an IBM-360 system was used to give some indication of a fast serial 
computation of the coefficients for various numbers of sample. The 
average time to perform an arithmetic operation (complex addition 
followed by complex multiplication) was estimated at 22 usee. From 
that average operation time, the maximum computation times for several 
different numbers of samples were computed and are shown in Fig. 6. 

The question of combining the speed of parallel processing with 
the simplicity of serial operation must be answered with a compromise. 
It is quite evident that for 16 samples or more, straight parallel 
computation is unreasonably expensive. With a very low hypothetical 
price of $10,00 per adder, a 16 sample processor’s arithmetic unit 
alone would cost $20,000.00. Consider, however, a unit which would 
compute 16 coefficients according to the flow in Fig. 4 by using a 
single full-parallel eight complex sample processor and serially 
processing the basic couplet’s arithmetic and complex multiplications. 

A general 12 x 12 bit CSA multiplier that will give the product 
in 360 nsec may be constructed with 132 adders and 156 AND gates. 

With four such multipliers and two simple adders a complex multiplier 
is constructed that will provide the product of two complex quantities 
in 380 nsec. Since the basic couplets characteristically fall into 
sum and difference couplets (see the first pass column in Fig. 2) it 
seems reasonable to compute both at the same time. By adding a complex 
subtractor (two words) before the multiplier and a parallel adder, the 
sum couplet and the difference couplet with complex multiplication are 
computed in 400 nsec. Such a unit would require approximately 575 



29 



adders and 624 AND gates. If 9 AND gates are considered equivalent 
to one full adder, the above circuit is equivalent to 640 adders. 

At this point the *'hybrid” processor consist of one eight- 
comp lex- sample processor and one arithmetic unit that computes 
the sum couplet, the difference couplet, and multiplies the dif- 
ference couplet by a complex constant in a single 400 nsec operation. 

For the 16 sample problem of Fig. 4, the arithmetic unit would 
only be required to cycle through eight calculations (First Pass), 
then the eight sample unit would cycle twice. Assuming a .25 
usec read and write time, one complete arithmetic operation could be 
completed in less than 1.25 usec. The eight arithmetic operations 
would require 10 usec and the two cycles of the eight sample unit 

9 usec (4.5 usec per cycle) for a 19 usec computation time. The 

above calculation times have been extended for various numbers of 
sampled inputs and the results are shown in Fig. 6. 

The hybrid processor *s arithmetic unit at this point consist of 
a 1240 adder unit. However, the processor must have a control unit 
and a control memory. The control memory requires a maximum of N/2 
complex multipliers plus the control program. The control unit must 
be able to index the pass numbers, the arithmetic operations, and the 
eight sample process pass number. Assuming a single word in memory 
is equivalent to a full adder in cost and complexity, a possible basis 

of comparison is achieved. Since two words of memory are needed for 

each multiplier there will be N words or, equivalently, N adders 
required in the control memory of a N sample unit. An assumption of 
a 200 adder equivalent cTontrol unit, including the control program 
memory is a reasonable estimate. For the range of sample inputs from 
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8 to 1024 this value should remain relatively constant. In general, 
the control unit then requires the equivalent of N + 200 adders. 

The 16 sample hybrid processor is the equivalent of 1450 full 
adders. Figure 7 illustrates the relative complexities of the hybrid, 
the parallel, and the serial processor. The basic requirements for 
sample input and coefficient output memory is the same for all three 
processes and is not included in the complexity figures. 
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FIG. 6. COMPUTATION TIME Vs. NUMBER OF SAMPLES 
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COMPLEXITY (Full-Adder Equivalents) 




FIGURE 7. COMPLEXITY Vs. NUMBER OF SAMPLES 
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4. Conclusions 



Figure 6 graphically shows the superiority of the parallel 
processor with respect to speed of computation while Fig. 7 demon- 
strates the enormity of size and complexity required to obtain those 
speeds. In reverse, Fig. 7 shows the serial processor with the 
complexity advantage and being extremely slow in computation, as 
shown in Fig. 6. 

The combining of serial and parallel operations succeeded in 
incorporating the better qualities of each. Certainly, the com- 
plexity of the hybrid system follows the same trend as the serial 
processor in that it is relatively flat in the log-log plot. This 
was to be expected. In presenting the hybrid system, the major 
nature of the system is serial. The parallel characteristics are 
found in the handling of several computations at one time and reducing 
substantially the arithmetic operation time by batch processing 
through the eight sample full- parallel processor in the final stages. 

On the log- log plot of Fig. 6 the slope of the hybrid curve is slightly 
greater than that of the serial curve. However, at the last point 
presented (1024 samples) the hybrid processor is still 500 times 
faster than the serial unit. 

It is evident from the computation times shown that definite 
applications for the special purpose handling of FFT*s do exist. In 
fact, for many applications, real-time analysis is possible. A central 
problem remains, that of analog to digital conversion of the data in 
times that will be able to fully utilize the speed of the processors. 
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APPENDIX I 



EQUATIONS FOR EIGHT SAMPLE APPLICATION 
The following set of equations defines the DFT for n=8. The 
equations are expressed in terms of the sample number (sample subscript) 
rather than using the subscripted notation of the Introduction. In 
this notation the couplet [Z(0)+Z(4)] would appear as (0+4), and non- 
integer numbers are multiplicative constants. 

F(0) = [(0+4)+(2+6)] + [(l+5)+(3+7)] 

F(l) = U0-4)+.707[(l-5)-(3-7)] j -j j (2-6)+.707[ (l-5)+(3-7)] 1 
F(2) = [ (0+4)- (2+6)]- j[ (1+5)- (3+7)] 

F(3) = j(0-4)-.707[(l-5)-(3-7)] J +j j (2-6)- . 707[(91-5)+(3-7 ) ] ’ 

F(4) = [ (0+4)+(2+6)]-[ (l+5)+(3+7)] 

F(5) = (0-4)-.707[(l-5)-(3-7)l j -j j (2-6)- .707[ (l-5)+(3-7)] , 

F(6) = [(0+4)-(2+6)]+j[(l+5)-(3+7)] 

F(7) = I (0-4)+.707[ (l-5)-(3-7)] j +j i,(2-6)+.707[ (l-5)+(3-7)] j 
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APPENDIX II 



EQUATIONS FOR SIXTEEN SAMPLE APPLICATION 
The set of equations beginning on the following page defines 
the DFT for N=16. As in Appendix, I, the equations are expressed 
in terms of the sample number. 



37 



F(0) = i[(0+8)+(4+12)]4{ (2+10)+(6+14)] +i [ (l+9)+(5+13)]4{ (3+ll)+(7+15) ] 

i 

F(l) = [ i (0-8)+.707[(2-10)-(6-14)] j +.92 ( 1- 9)+.707[ (3- 11)- (7- 15)] , -.38.* (5- 13)+.707[ (3- ll)+(7- 15)] j 

-j[[ (4-12)+.707[(2-10)+(6-14)] j +.38 j (l-9)+.707[ (3- 11)- (7- 15)] | +.92 j (5-13)+.707[(3-ll)+(7-15)] | 
F(2) = [[(048)- (4+12)] + . 707 j [ (1+9)- (5+13)]-[ (3+11)- (7+15)] 

-jf[(2+10)-(6+14)+.707 j [(l+9)-(5+13)]+[(3+ll)-(7+15)]' 5 
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F(10) = ■*[ (0+8)- (4+12)]-. 707 I [(1+9)- (5+13)]-[ (3+11)- (7+15)] 1 
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APPENDIX III 



ILLUSTRATIVE EXAMPLE OF THE DISCRETE 
FOURIER TRANSFORM AND ITS INVERSE 
A time signal, ^(t), is sampled eight times at equally spaced 
intervals. For convenience, the signal is a sinewave and the 
interval tt/4 is chosen. 




The sampled signal appears as 



Z(k) 

1 - 
.707 - 



! ° 

I 

-.707 -! 

-1 - 




6 7k 



Z(0)=Z(2)=.707 Z(l)=l 

Z(4)=Z(6)=- ,707 Z(5)=-l 

Z(3)=Z(7)=0 



Forming the basic couplets of the equations in Appendix I. 



Z(0)+Z(4)=0 

Z(2)+Z(6)=0 

Z(l)+Z(5)=0 

Z(3)+Z(7)=0 



Z(0)-Z(4)=1.414 

Z(2)-Z(6)=1.414 

Z(l)-Z(5)=2 

Z(3)-Z(7)=0 
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Now, solving for the coefficients 
F(0)=(0+0)+(0+0) = 0 

F(l)=[l.414+.707(2-0)]-j[l.414+. 707(2+0)] = 2.828(l-j) 
F(2)=(0-0)-j(0-0) = 0 

F(3)=[1.414-. 707(2-0) ]+j[l.414-. 707(2+0)] = O+jO = 0 
F(4)=(0+0)-(0+0) = 0 

F(5)=[1.414-. 707(2-0)]- j[l.4l4-. 707(2+0)] = O-jO = 0 
F(6)=(0-0)+j(0+0) = 0 

F(7)=[l.41+. 707(2-0) ]+j[l.414+. 707(2+0)] = 2.828(l+j) 

If the foregoing procedure is correct, application of the inverse 
transform should produce the original sample values. Noting that the only 
non-zero coefficients are F(l) and F(7), the inverse is written 
Z (k)=l/8 F(l)exp( jnk/4)H-F(7)exp( jnk/4) 
and the Z(k)'s are 

Z(0)=(2.828/8)[(l-j)+(l+j)] = .707 

Z(l)=(2.828/8)[(l-j)(.707)(l+j)+(l+j)(.707)(l-j)] = 1 
Z(2)=(2.828/8)[(l-j)(j)+(l+i)(-j)] = .707 
Z(3) = (2.828/8)[(-.707)(l- j)^+(- .707) (1+j)^] = 0 
Z(4)=(2.828/8)[(l-j)(-l)+(l+j)(-D] = -.707 
Z(5)=(2.828/8)[(-.707)(l+j)(l-j)+(-.707)(l-j)(l+j)] = -1 
Z(6)=(2.828/8)[(l-j)(-j)+(l+j)(j)] = ^.707 
Z(7)=(2.828/8)[(.707)(l-j)^+(.707)(l+j)^] = 0 

The above Z(k) match precisely the original data; thus showing the 
correctness of the development 
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APPENDIX IV 



CARRY- SAVE ADDERS 

Basically a full adder has three inputs (two addends and a 
carry-in) and two outputs (sum and carry-out). Thus, it may be 
looked upon as a reduction process of from three to two arguments. 
Considering the "carry- in" as a third addend, it is apparent that 
the addition of more than two quantities may be handled by cascading 
the afore mentioned reduction process. Such implementation of full 
adders has led to both the process and the adders themselves being 
called Carry-Save Adders (CSA) . 

If more than three addends are being summed more than one 
initial CSA would be required. The figure below schematically shows 
the reduction process for a six argument addition. 




Schematic representation of a six argument Carry-Save Addition 

Note that this is strictly a representation of the reduction process. 
In the actual process the "carry-out" of each adder must go to the 
next higher significant bit train. Its position, as shown in the 
diagram, would be taken by the carries from the next lower significant 
bit train. It is noted here that the CSA process generally requires 
k-1 units for the reduction of k addends. The final stage (level) of 
the CSA process is of interest since the propagation of the carries 
are of a "ripple" nature. 
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The following diagram presents a complete four-bit adder 
for four addends (A,B,C, & D) . It is obvious that the CSA process 



0 

V 

E 

R 

F 

L 

0 

W 




Block diagram of a four-bit, four addend, 
Carry-Save Adder (an addition of four 4-bit words) 



is strictly a logic function. With the assumption of a realistic 
signal delay time through each logic level, the total computation 
time may be calculated. Although pressing the state of the art 
at this time, an assumption of a 20 nsec full adder is not unrealistic. 
Thus, for the given adder, three levels of delay plus the final carries* 
ripple through two delay levels gives a total of five delays, or a 100 
nsec addition of four, 4-bit arguments. 

With increasing word lengths the computation time increases 
linearly for the CSA process. This characteristic is the result of 
adding One more stage of carry ripple for each additional bit. For 
a 12-bit version of the four argument adder, there would be 13 units 
of delay and an addition time of 260 nsec. 
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The use of carry look ahead techniques can reduce the ripple 
time by a factor of two per stage of delay. For the 12-bit example, 
the three initial units of delay are summed with 10 half-units for 
8 delays and a 160 nsec addition of four, 12-bit, quantities is 
possible . 

Finally, ones complement arithmetic is appropriate since it is 
easily obtained by using the two's complement and just adding one 
in the least significant bit train. 

The following table gives the number of delays as a function of 
the number of arguments in a CSA process. 



Number of Arguments 

1 ... . 

3 . . . . 

4 . . . . 

5,6 . . . . 

7-9 . . . . 



10-13 

14-19 

20-28 



Number of Delays 

... 1 
... 2 
... 3 

... 4 

... 5 

... 6 

... 7 

... 8 
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